Entanglement spectroscopy of a driven solid-state qubit and its detector 
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We study the asymptotic dynamics of a driven quantum two level system coupled via a quantum 
detector to the environment. We find multi-photon resonances which are due to the entanglement 
of the qubit and the detector. Different regimes are studied by employing a perturbative Floquet- 
Born-Markov approach for the qubit+detector system, as well as non-perturbative real-time path 
integral schemes for the driven spin-boson system. We find analytical results for the resonances, 
including the red and the blue sidebands. They agree well with those of exact ab-initio calculations. 
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A prominent physical model to study dissipative and 
decoherence effects in quantum mechanics is the spin- 
boson model [jj. Currently, we witness its revival since 
it allows a quantitative description of solid-state quantum 
bits (qubits) . A more realistic description requires the 
inclusion of the external control fields as well as the de- 
tector. In the spin-boson model, the environment is char- 
acterized by a spectral density J(u). In its widest used 
form, J(lu) is proportional to the frequency ui mimick- 
ing the effects of an Ohmic electromagnetic environment. 
However, if the environment is formed by a quantum de- 
tector which itself is damped by Ohmic fluctuations, the 
simple Ohmic description might become inappropriate. 
As an example, we focus on a superconducting ring with 
three Josephson junctions (so termed flux-qubit). It is 
read out by a dc- SQUID 0- whose plasma reso- 
nance at Q p gives rise to an effective spectral density 
J G ff(cj) for the qubit with a peak at Q, p 6], cf. Eq. |@J be- 
low. Recently, the coherent coupling of a single photon 
mode and a superconducting charge qubit has also been 
studied Q. Until now, the effects of such a structured 
spectral density on decoherence and in presence of a res- 
onant control field have only been studied in 0, [|J within 
a perturbative approach in J e g. It was shown in (lfjflU 
for the static case that a perturbative approach breaks 
down for strong qubit-detector coupling, and when the 
qubit and detector frequencies are comparable. 

In the presence of microwaves, multi-photon reso- 
nances are expected to occur when the frequency of the 
ac-ficld, or integer multiples of it, match characteristic 
energy scales of the system [l2j. Such multiphoton res- 
onances can be experimentally detected in an ac-driven 
flux qubit by measuring the asymptotic occupationprob- 
abilities of the qubit, as the dc-field is varied ll.3| . 
These "conventional" resonances, which have also been 
theoretically investigated in |14). could be explained in 
terms of intrinsic transitions in a driven spin-boson sys- 
tem with an unstructured environment. 

In this Letter, we investigate the asymptotic dynamics 
of a quantum two state system (TSS) with a structured 
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FIG. 1: Schematic picture of the models we use. In (a) the 
system is a two-level-system (TSS) coupled to a harmonic 
oscillator with the latter coupled to an Ohmic environment 
with spectral density Johm(w). In (&), the TSS is coupled to 
an environment with peaked spectral density J e ff (w). 



environment, simultaneously driven by dc- and ac-fields. 
We show that a strong coupling between qubit and detec- 
tor, together with the presence of a control field, yields a 
non trivial dynamics involving additional resonances in 
the entangled qubit+detector system. Our results are in 
agreement with recent experimental findings where such 
"unconventional" multi-photon transitions have been ob- 
served 0. We evaluate the TSS dynamics in two com- 
pletely equivalent models, cf. Fig. 1. In model (a), the 
TSS is coupled to its detector being represented as a sin- 
gle harmonic oscillator (HO) mode with frequency £l p 
with interaction strength g. The HO itself interacts with 
a set of harmonic oscillators, cf. Fig. la. The correspond- 
ing Hamiltonian is -ffQos(^) — HQo{t) + Hob-, where 
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Here, <Ji are Pauli matrices, fiA is the tunnel splitting, 
and e(t) = e + scos(ftt) describes the time-dependent 
driving with the static bias £o- For s = 0, the level split- 
ting of the isolated TSS is hv = HWEq + A 2 . Moreover, 
B is the annihilation operator of the localized HO mode, 
X = B* + B, while b k denote the bath mode opera- 
tors. The spectral density of the continuous bath modes 
is Ohmic with dimcnsionlcss damping strength k, i.e., 



^Ohm(w) = ^ " = KUJ ~ 



(2) 



where we have introduced a high-frequency cut-off at u>d- 
In this approach, we shall consider the combined TSS + 
HO as the central quantum system. 

In the second model (b), we exploit the exact one-to- 
one mapping yjj| of the Hamiltonian (JJJ onto that of a 
driven spin-boson Hamiltonian 0] 
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where b k is the annihilation operator of the k— th bath 
mode with frequency £j k ■ Following |6( , the spectral den- 
sity has a Lorentzian peak of width 7 = 2irnQ p at the 
characteristic detector frequency f2 p . It behaves Ohmi- 
cally at low frequencies with the dimensionless coupling 
strength a = lim w _>o J e s(u>) /2u> and reads 



u> 



2\2 



(4) 



The relation between g and a follows as g — £l p \J a/8n. 
In this model, we associate the detector as part of the 
qubit environment. 

The qubit dynamics is described by the reduced den- 
sity operator p(t) obtained by tracing out all environ- 
mental degrees of freedom. We study the population dif- 
ference P(t) := (a z )(t) = tr(p(t)a z ) in the asymptotic 
limit, i.e., P x = lim t _» 00 (P(i))n, where the averaging is 
over one period of the ac-field. 

Case of weak damping and low temperatures. For 
k < 1 and ksT < SA, it is convenient to use model 
(a). The equations of motion for the TSS+HO re- 
duced density matrix are most conveniently derived in 
the Floquet basis The Floquet states \<j> a (t)) — 

Y] \4>^) exp(infii) corresponding to a periodic Hamil- 
tonian H(t) can be obtained from the eigenvalue equa- 
tion H\(j> a (t)) = s a \<j> a (t)), with the Floquet Hamiltonian 
TL = H(t) — ih-^. Upon including dissipative effects to 
lowest order in n, a Floquet-Born-Markov master equa- 
tion is obtained [l^lljj. We average the 27r/f2-periodic 
coefficients over one period of the driving, assuming that 
dissipative effects are relevant on much larger timescales. 
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FIG. 2: Left: Quasi-energy spectrum e atk of the driven 
TSS+HO system vs dc-bias £0 (in units of A). The quasi- 
energies are defined up to an integer multiple of hQ, i.e., 
£ a ,k = £ a + khfl. Inset: Zoom of an anti-crossing. Right: Poc 
exhibits sidebands corresponding to quasi-energy level anti- 
crossings. Parameters are f2 = 10A, s = 4A,g = 0.4A, Q p = 
4A, k = 0.014 and k B T = O.lftA. 



In the Floquet basis, this yields equations of motions for 
p a p(t) = ((j) a (t)\p(t)\cl>p(t)) of the form 



Qap(t) = ~j( £a ~ e fi)Qap{t) + ^2 L a p,a'P'Qa>p>(t), (5) 

a'f3' 



with the dissipative transition rates 

L a p,ct'P> = y^(-^aa',n + Nppi , n )X aa > , n Xp' p,-n 
n 

— Spp 1 2_j Np"a' ,nX a p" ,-nXp» a ' ,n 

P",n 

— S aa > N a iip> tn Xpi a n t - n X a 'ip >n . (6) 



^\X\4,f +n) ), and N aP , n = N(e a 
2h \ 



(coth ( 2k £ BT ) — I) (assuming 



Here X a p^ n = Y,kW< 
E/3 + nW.) with N(e 

LU£, — > OO). 

Following Ref. 01 we write the Floquet Hamiltonian 
TLqo m the basis |o, n), with \a) = \g/e,m), g/e being 
the ground/ excited state of the qubit, to the oscillator 
state, and n the Fourier index. In this basis, TLqo has 
diagonal elements TL a n,an = h[=fi//2 + mQ p + nfi], and 
off-diagonal elements V an M = (a\S njk hgXo- z + (5 n ,k+i + 
&n+i,k)^fo- z \b) . The quasi-energy spectrum of TLqo is 
shown in Fig. 2 as a function of the bias eo- We 
find avoided level crossings when E an ^ m := TL an ,an — 
TLbmfim = + 0(V 2 ), i.e., (to > 0, -co < n< +oo) 

v = nn±mn p + 0(V 2 ), nQ = mflp + 0(V 2 ) . (7) 

Associated to the avoided crossings are resonant 
peaks/dips of Poo, see Fig. 2. The resonances at v = 
nfl ± mflp are known as red/blue sidebands [l9j |. 

In the following, we derive an analytical expression 
for the first blue sideband at v ps — fl p . Other res- 
onances can be evaluated in the same way. We include 
only one HO level (to = 0,1) which is appropriate be- 
cause we investigate a resonance between \g, 0) and |e, 1) 
with g/^lp <C f ■ We consider V an ,bk as a perturbation, 



3 



and use the method of Ref. 0, |2jj to obtain an effective 
Hamiltonian H c s = e lS H.Qoe~ lS , with 
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for \E an . bm \ ^ \v + Slp-Sl\, (8) 



and iS a 



'an,bm — for \Ean,bm\ |'' ' --i> 

block-diagonal H c s has the same eigenvalues as TIqo 
with quasi-degenerate eigenvalues £1,2 in one block. With 

C V3 = SC^T^Tgp^&F and <5 = «/-n+fi f ,-2c 1 
the quasi-energies up to second order in V read 



ni. The 



ei/a = -W2 + M/2(1t Vl + Ai/<5 2 ) + fici, 
£3 = —Hi'/2 + hilp — Hc3, £4 = ft^/2 + ftc3. (9) 

The eigenvalue splitting at the level crossing, <5 = 0, is 

hAe gs [n 2 + n 2 p + + n p )] 



4i/(n - i/)nn p (i/ + Q, r 



(10) 



The Floquet states are, with tan# = 2|Ai|/<5, B + (x) 
cos(a;) and B~{x) = sin(x), 
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-iS 



5,1), 







-iSi 



!,0). 



(11) 



With this, we can calculate the rates in Eq. © up to 
second order in V. It holds, L1122 = £2211 = £3344 = 
£4433 = 0. To find the stationary state of Eq. (0), we 
assume that p a p (00) = for a ^ j3, except for px2 and 
P21 (secular approximation). This is valid if e a — ep 3> 
L a p. a 'f3', which is true for non-quasi-degenerate eigenval- 
ues because k <C 1. We find at resonance 



(12) 



implying a complete inversion of population at low tem- 
peratures. Far enough off-resonance, we can assume that 
Pi2(oo) = £21(00) = and sin (9/2) ~ 9/2. We pre- 
sume ksT <C Ml p , Ml, hv (which allows to set N(Ml p ) — 
N(MY) = N{hv) — 0) and find the major result 
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|.(13) 



Because the oscillator can give its energy directly to the 
environment, the decay from \e/g,l) to \e/g, 0) is much 
faster than the other processes and does not play a role 
in Eq. (|13|) . Hence, P^ is determined by the ratio of 
two rates: L4411 ~ sin 2 (9) ~ s 2 g 2 which describes the 
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FIG. 3: First blue sideband Poo vs eo (in units of A) at 
v = Q — fi p . The solid lines are the analytical prediction 
O for (a) g = 0.05A, (6) g = 0.2A, (c) 3 = 0.4A. The tri- 
angles are the results of a Floquet-Bloch-Redfield simulation, 
cf. Eq. ijSy. with one (upward triangles) and two (downward 
triangles) HO levels taken into account. The circles in (b) are 
the results from a QUAPI simulation with six HO levels. We 
choose s = 2A, fi = 10A, n p = 4A, « = 0.014, fc s T = O.lfiA. 



timescale of driving induced transitions from \g, 0) to 
|e, 0), and Z/1144 ~ g 2 for the qubit decay from |e,0) 
to \g, 0) via the oscillator. Since both scale as g 2 we 
find for this particular resonance that is indepen- 
dent of g. Fig. |3| shows the results of Eq. I|13[) and dif- 
ferent numerical results, including those of an ab-initio 
real-time QUAPI [2l| calculation. A good agreement, 
even near resonance, is found. A similar analysis yields 
Poo = f tanh (^r) + 0(V 2 ) for the first red sideband 
at v = fl + Qp, which is very close to thermal equilibrium 



for low T. For v = Q only the oscillator is excited and 
thermal equilibrium P^ = ^ tanh( 2 ^ T ) is recovered. 

Case of strong damping and/or high temperatures. In 
the complementary regime of large environmental cou- 
pling and/or high temperatures it is convenient to employ 
model (6), and is appropriate to treat the system dynam- 
ics within the noninteracting-blip approximation (NIBA) 
0. The NIBA is non-perturbative in the coupling a but 
perturbative in the tunneling splitting A. Within the 
NIBA, and for large driving frequencies £1 > A, one finds 
-Poo = K l k a .12], where 



dTh ± (T)B ± (e T)J Q 
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(14) 



The influence of the dc- and ac-field is in the terms 
B ± (x), and in the Bessel function J , respectively. Dissi- 
pative effects are captured by h ± {t) = e- Q '^B ± [Q"(t)}, 
where Q'(t) and Q"(t) are the real and imaginary parts 
of the bath correlation functional . For the peaked spec- 
tral density Eq. (0J one finds 

Q'{t) = Q' 1 (i)-e- rt [Y 1 cos(n p <) + T 2 sin(n p i)] 
Q"{i) = A x - e~ rt L4i cos(O p i) + A 2 sin(f2 p t)] . (15) 

Here, (3 = h/k B T, T = irnflp, Cl 2 p = fi| - T 2 and 
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FIG. 4: Poo vs £ (in units of A). The solid line is the NIBA 
prediction, while the circles are from a QUAPI simulation 
with 6 HO levels (g = 3A, s = 4A, Q = 10A, k = 0.014, 
k B T = 0.5ftA, tt p = 4A). Inset (a): NIBA result for k B T = 
2hA. The arrows indicates the first red sideband at v = 
O + Q p . Inset (b): Q'(t) vs t shows damped oscillations. 



we could not reach the parameter regime T/fl p <C 1, 
but still clear resonance dips are observed at eo = ^) 
Eq = — flp and eg — Q — 2tt p . For fc^T ~ M7 p , we also 
find the first red sideband a,t e$ = £l + £l p , see inset. 

In conclusion we evaluated the asymptotic population 
of a driven TSS in a structured environment. We have 
derived analytic expressions for the shape of the reso- 
nances for both weak and strong damping. We show that 
the coupling of the TSS and the detector is revealed in 
the occurrence of characteristic multi-photon resonances, 
also reported in recent experiments |5j , in the asymptotic 
population of the TSS. A complete population inversion 
is predicted for the blue-sidebands transitions, while val- 
ues close to equilibrium are found for the red-sidebands. 

Support by the Dutch NWO/FOM, and the Univer- 
sitatsstiftung Hans Vielberth, and discussions with P. 
Bcrtet, I. Chiorescu and H. Mooij are acknowledged. 



where v n = 2nn/f3. Moreover, C — cosh(/9f2 p ) — cos(/3r), 
CY 1/2 = tA 2 /i sinh pCl p - A l/2 sin f3T, A 2 = aTr(T 2 - 
Cl p )/2T0, p , Ai = na. So, Q' and Q" display damped os- 
cillations (cf. Fig. not present for a pure Ohmic spec- 
trum. It is the interplay between these oscillations and 
the driving field which induces the extra resonances in 
Poo. In the regime T/il p <C 1, the term exp(— Tt) in Eq. 
(|16fl varies slowly on the time-scale of the oscillations. 
We expand Q 1 and Q" as well as the Bessel function Jo 
entering Ijl4(l using Bessel function identities and find the 
important result 

k o = E E W dte-^f± n (t) , (17) 

Tci— n— — oo 

where e mn = eo — mCl p — ?if2, and 
Re 

fmn(t) = t [c™„(i)cos(e mn i) ± c* n {t) sin(e m „t)] , 

4« = Jn (~) J m (e- rt c 1 ) J B ± (m0)(-*) m e-^ 1 . (18) 

Here is w\ = x/C^i - iYi) 2 + (A 2 - iY 2 ) 2 , and tan tj> = 
-(A 2 -iY 2 )/(A 1 -iY 1 ). Thus, from Eq. (EJ we ex- 
pect resonances when e nrn = 0. Without driv- 
ing we always find that around Eq = mfl p it holds 
Poo w tanh(m/3i7p/2), since limr/n — >otan(m^>) = 
itanh(m/3f2 p /2) (for not too large T, i.e., cos(/?r) <C 
cosh(/351 p )). Hence, Poo acquires its NIBA thermal 
equilibrium value, and driving is needed to see reso- 
nances. For "conventional" resonances at £o — n ^ w e 
find Pop ~ 0, as predicted for unstructured environ- 
ments [rj, |22j. Finally, for eo = nfl ± m£l p , we recover 
Poo ~ ± tanh(m/3f2 p /2), as also was found within the 
Floquet-Born-Markov approach, cf. (|12(l . Results of a 
numerical evaluation of Poo are shown in Fig. 4, using 
the NIBA result (|18fl . as well as the exact ab-initio real- 
time QUAPI method |2l|. In the numerical evaluation, 
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